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In columnar assemblies of helical bio-molecules the azimuthal degrees of freedom, i.e. rotations 
about the long axes of molecules, may be important in determining the structure of the assemblies 
' especially when the interaction energy between neighbouring molecules explicitly depends on their 

' relative azimuthal orientations. For DNA this leads to a rich variety of mesophases for columnar 

, assemblies, each categorized by a specific azimuthal ordering. In a preceding paper [A.Wynveen, 

■ D. J. Lee, and A. A. Kornyshev, Eur. Phys. J. E, 16, 303 (2005)] a statistical mechanical theory 

^ ' was developed for the assemblies of torsionally rigid molecues in order to determine how thermal 

^ , fluctuations influence the structure of these mesophases. Here we extend this theory by including 

torsional fluctuations of the molecules, where a DNA molecule may twist about its long axis at the 
cost of torsional elastic energy. Comparing this with the previous study, we flnd that inclusion of 
torsional fluctuations further increases the density at which the transition between the hexagonal 
structure and the predicted rhombic phase occurs and reduces the level of distortion in the rhombic 
phase. As X-ray diffraction may probe the 2-D lattice structure of such assemblies and provide 
information concerning the underlying interaction between molecules, we have also calculated cor- 
relation functions for the azimuthal ordering which are manifest in an x-ray scattering intensity 
^ ' profiles. 

O 

' I. INTRODUCTION 

' Condensed DNA exists in a rich variety of phases and mesophases. Molecular assemblies of DNA are considered to 
be lyotropic, so much of this rich phase structure may be accessed by changing the concentration of DNA in solution. 
As the concentration of DNA molecules is increased, the orientations of the molecules evolve from a completely 
disordered isotropic phase into a liquid crystal phase. Eventually at large enough densities, long range order is 
^ ■ established, and a crystalline state is established jj. In the liquid crystal phase, various mesophases are observed 
0, 0, 0, ■ At relatively low concentrations the DNA molecules form a cholesteric mesophase, whereas at larger 
concentrations the molecules form columnar assemblies. We should point out that not only is DNA concentration a 
parameter in determining what mesophase the DNA is in, but other factors such as monovalent salt concentration, 
T-H ' as well as the type and quantity of condensing agent used, are equally important Hence, any analytical study 

' is necessarily quite involved due to the complexity system. 

O I 

c/j ' Many of these mesophases have been seen in biological systems |^ . Determining how these mesophases form is 
O , important for understanding how DNA packs into viral capsids and sperm heads jQ] and may also be relevant in 
c/2 ' the advancement of gene therapies 0,01 ■ The study of x-ray diffraction patterns of such mesophases, as well as 
an understanding of the statistical physics underlying DNA assemblies, could well elucidate the nature of the forces 
• between such molecules llj. Furthermore, understanding how properties of these mesophases might depend on the 
sequence of base pair text of the DNA that form the mesophase may reveal information concerning the process of 
homologous recognition of genes p^ . 
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Another question that might be probed in the study of DNA mesophases is whether the interaction between 
' molecules depends on the azimuthal orientations of the DNA about their long axes. The polyelectrolyte model 
[r^lT3 . ll5j assumes that the forces between molecules are independent of the azimuthal orientations of the molecules. 
This is only valid if the azimuthal dependence of the force is very weak or is completely screened by the solvent and 
its constituent ions. This model, however, is insufficient to describe the various mesophases. An alternative model 
[T^ was proposed, which took into account the helical nature of the surface charge patterns on the DNA molecule. 
In this model the interaction energy of a pair of molecules depends strongly on the azimuthal orientation out to 
separations of i? > 30A. This azimuthal dependence was further enhanced when counterions were assumed to be 
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preferentially adsorbed within the grooves of the molecule |17|. Such an azimuthal dependence was indeed shown to 
influence both the structures of columnar and cholesteric assemblies Il9l. 



In Rcf. 18] it was shown that the variety of states observed in columnar assemblies are intrinsically linked 
with the azimuthal orientation of the molecules. These states corresponded to different "spin"-orderings, i.e. the 
configuration of the relative azimuthal orientations between neighboring molecules, of the molecules situated on the 
two-dimensional lattice of a columnar assembly. For example, there exist "ferromagnetic" and "antiferromagnetic" 
states, i.e. states corresponding to where all the molecules are azimuthally aligned and to where molecules along 
a lattice direction have alternating values for their azimuthal orientations. In a later work we extended this 
ground state calculation |0 to incorporate the effects of thermal fluctuations to build up a full statistical mechanical 
model of the columnar assemblies. Here, two new transitions were observed for a fixed hexagonal lattice. The 
first corresponded to a transition from one of two topologically distinct states of a three-spin configuration, where 
nearest neighbors of a molecule have one of two different azimuthal orientations (referred to as the Potts state), to 
a more disordered state where both topologies were equally likely in the system. The second was a transition to a 
Berezinskii-Kosterlitz-Thouless like vortex state '2l|. Here, also, it was shown that the antiferromagnetic state was 
only stable when the two-dimensional hexagonal lattice was distorted. Allowing for these lattice distortions, a phase 
transition occurring between this rhombic (distorted hexagonal) antiferromagnetic state and the ferromagnetic state 
is observed. As compared to the ground state calculations, incorporating thermal fluctuations resulted in a decrease 
in the mean separation between molecules at which this transition occurred, a reduction in the amount the distortion 
to the hexagonal lattice, and a shift from a second order transition to a first order one. In both of these studies 
[T8Ll20l| . however, the DNA molecules were assumed to be completely rigid. 

Torsional flexibility has been considered in a ground state calculation of interacting nonhomolgous DNA 
where the reduction of sequence-dependent distortions due to electrostatic interactions between the molecules is 
observed. This present study, however, considers the effects of thermal excitations of the flexible molecules so 
that torsional fluctuations, in which the azimuthal orientation of the double helix is no longer uniform along 
the length of the molecule, occur. Such effects may be implemented relatively easily into an energy functional, 
which forms a starting point for the statistical mechanical treatment. Previously, for rigid molecules, the effect of 
doubling the length has the same effect as halving the temperature. In flexible molecules, however, length is no 
longer such a trivial parameter in the theory. Here, shortening the molecules plays a similar role to increasing the 
torsional rigidity since the free energy cost of twisting a short molecule to the same degree, i.e. the same variation 
in the azimuthal angle between the ends of the molecule, of that of twisting a longer molecule is much greater. 
Hence, very short molecules may be assumed to be rigid. Furthermore, as a molecule is lengthened beyond a certain 
torsional persistence length, the extent of the torsional fluctuations becomes independent of the length of the molecule. 

Though formulating the energy functional and partition function for columnar assemblies of flexible molecules is 
a relatively simple task, calculating the free energy and other thermodynamic quantities is not. Due to the extra 
degrees of freedom associated with torsional fluctuations, Monte-Carlo methods become less reliable and more 
time consuming. Utilizing self-consistent approximations developed through field theoretical methods, however, we 
encounter a problem arising from treating molecules of finite length; how to account for the freely fluctuating ends of 
the molecules in the assembly. Nevertheless, we have managed to incorporate free end effects into the field theoretical 
calculations and have determined the extent to which they alter the results. 

The layout of the paper is as follows. In Sec. II we consider the case of finite temperature interactions between 
two molecules in parallel juxtaposition as a starting point for developing the present calculations since the reliability 
of the these calculations may be tested against the "Quantum Mechanical" formulation of DNA interactions set 
out in Ref. |2^. In Sec. HI we summarize the results of the various ground state calculations of the system and 
then incorporate thermal fluctuations at the harmonic level (Gaussian fluctuations) for the different configurations 
of the columnar assembly. Here, we determine the free energy of the assembly as well as the correlation function 
corresponding to the variation of relative azimuthal fluctuations between molecules which is important in determining 
the x-ray diffraction patterns of the assemblies IflJ. In Sec. IV we go beyond the harmonic treatment of the 
interaction, developing a self-consistent approximation for the assembly. This relies on a Hartree approximation for 
the case where the ends satisfy periodic boundary conditions, i.e. both ends of a molecule have the same azimuthal 
orientation, and a correction which takes into account independent fluctuations at the ends of each molecule. 
At this level of the calculation, we demonstrate how torsional fluctuations affect the free energy and correlation 
functions of the assemblies and compare these results to those of assemblies composed of rigid molecules |23|. We 
find that the correction accounting for independent fluctuations of the ends of the molecules is quite small, and so 
in first approximation may be neglected, demonstrating that the utilization of periodic boundary conditions in the 
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calculation is a reasonable approximation. Finally, in Sec. V, we discuss our results and possible future developments. 



II. HOW TO TREAT MOLECULES OF FINITE LENGTH: THE 2-BODY PROBLEM 

Our primary aim is to develop a field theoretical framework in which to treat the statistical mechanics of flexible 
molecules condensed in columnar assemblies. As mentioned in the introduction, we first encounter the problem of 
how to treat molecules of finite length. If we assume that periodic boundary conditions can be applied to the ends of 
the molecules, however, we may take advantage of field theoretical techniques that already have been developed 24J. 
By assuming periodic boundary conditions we presuppose that both ends of a DNA molecule fluctuate torsionally 
in phase. In other words, both ends share the same azimuthal orientation, which isn't necessarily the case. And so 
within our statistical mechanical treatment, we also take into account freely fluctuating ends of the molecules. 

Before turning our attention to an assembly, we begin by first considering two molecules in parallel. One of the 
main motivations for doing this is that for the pair interaction, an alternative approach for formulating the statistical 
mechanics of interactions between flexible molecules, which resembles that of a quantum mechanical problem, may 
be employed This, then, provides a check on the reliability of the formulation with which we use to treat the 
assemblies. 

The pair potential is a function of the interaxial separation R and the relative azimuthal angle (f> between molecules. 
The relative azimuthal angle is defined as the difference in the azimuthal angles between the two molecules, = (f>i — (f>2, 
where is the angle that a vector, from the center of the ith molecule to the middle of its minor groove, makes with 
an axis that lies perpendicular to the long axes of the molecules. This axis is chosen to pass through the centers of 
both molecules. The pair potential energy per unit length has the form 

Eir^tiR, </)) = ^ (-l)"a„(i?) cos(n0). (2.1) 

n=0 

Equation (2.1) is completely model independent and can be deduced purely from symmetry requirements. The 
first requirement is that rotation of one of the molecules one whole revolution about its long axis should leave the 
interaction energy unchanged. Hence we may express Eint as a Fourier expansion in terms of the relative azimuthal 
angle. Secondly, for DNA molecules, hehcal symmetry dictates that ii'int (-R, </>) = i?int(-R, — 0), so only cosine terms 
in the Fourier series are retained. Finally, provided that L ^ H, where H is the helical pitch of the DNA molecule, 
helical symmetry also ensures that the interaction energy per unit length does not depend on z, where z is the 
coordinate that runs along the long axis of the molecule. 

Before proceeding further, we need to consider the behavior of then a„ coefficients. We shall use forms for these 
coefficients obtained within the Kornyshev-Leikin (KL) theory of the electrostatic interactions between helical 
macromolecules. Such a treatment, grounded in Debye-Huckel theory, neglects important non-local effects p5ll2^ l27| 
of the dielectric response of explicit water. However, such a theory can still provide a qualitative picture into the 
behavior of our system. Furthermore, this theory should work quite well at sufficiently large enough separations 
between molecules and for dilute salt concentrations where Debye-Huckel theory is valid. Explicit forms for the a„ 
coefficients obtained in such a theory are given in Ref. ;23] . In the theory, all of the a„ coefficients decay exponentially 
with large interaxial separations with higher order terms (larger n) decaying at a greater rate. (The inverse decay 
lengths Kn for each these terms is given by Kn = \/ n^g^ -\- k^, where g = 2tt/ H and Ks is the inverse Debye screening 
length calculated in Donnan equilibrium [isl l2(l l2^.1 Because of this, it is sufficient to truncate the series given 
in Eq. (2.1) at n = 2 'l6]. The zeroth-order term oq, which describes the interaction between two cylinders with a 
uniform charge distribution, does not depend on (p and therefore need not be considered in the treatment of torsional 
fluctuations. However, it does play an important role in the positional structure of the assembly. Finally, since 02 
decays faster than ai, minimizing the interaction energy with respect to (p leads to a configuration below a critical 
value i?, of the interaxial spacing where the preferred value of (j) is non-zero, whereas above i?, it is zero. Even when 
including non-local effects, these qualitative features should still be manifest in the interactions. 

Upon including thermal fluctuations of the relative twisting of one DNA molecule with respect to another, (j) now 
must be assigned a z-dependence. Following from previous work |23i] . we write down the partition function as a 
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functional or path integral 

I?(/)(z) exp 



L/2 



dz 



-L/2 



dz 



ai cos((/)) + 02 cos(2(/)) 



(2.2) 



where C is defined as the torsional elasticity modulus of the helices. Here, we have assumed that the centers of both 
molecules are at z = and each molecule has length L. The free energy then can be found upon calculating the 
partition function. 

As was discussed previously in Ref. ^23'!, with this partition function the problem can be recast in an alternative 
formulation that mirrors a quantum mechanical problem. Here, we solve the Schrodinger equation 



where 



2 d(h^ 
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^('/') =-Tf cos(0) + ^cos(2<^), 



^0 



(2.3) 



(2.4) 



Ag — C/{2ai) and Ap — C/{2kBT). In this formulation, the free energy is expressed through the solutions of Eq. (2.3) 
as 



F = - fcsT In ^ (0+ ) exp ( - — j (0_ ) - fcsT In e 



EL\ 



(2.5) 



where is a constant which can be neglected when comparing free energies of different states. and 0_ are the 
values of (f) at the ends, z = L/2 and z = —L/2 respectively, of the molecules. Periodic boundary conditions amount 
to setting </)_|_ = (^_, and computing the free energy entails summing over However, as was pointed out before, 
both ends should be left to fluctuate freely. Consequently, we should integrate over both 0_ and 4>+, allowing (/)_ to 
take on any value between and 2tt, and 0+ to take on any value whatsoever. When the fluctuations are small we 
may expand V{(j}) to quadratic order in (j) around its preferred value in the ground state. We shall consider only the 
case when the preferred value is zero (large interaxial spacings), although this calculation may be easily extended to 
a non-zero preferred angle. Eq. (2.3) may be solved analytically to determine the free energy. Some of the details of 
the calculation of F are given in Appendix A. Here, we shall quote the end result for long molecules 



F 



-ksTQ +Eo 



CL 
4AA„ 



In 



keT f 2L\ 



(2.6) 



where Eq ~ —ai + 02 and A^ — C/{2{ai — 4a2)). The last two terms are leading order corrections due to the finite 
length of the molecules. 

Returning to Eq. (2.2), we now approach the problem directly through path integration, first in the Gaussian 
or harmonic approximation. In such an approximation, both cosine terms are expanded out to quadratic order in 
(j) around the preferred orientation. This is equivalent to expanding V^cjj) to quadratic order in (jj in the previous 
formulation. We have demonstrated this equivalence in Ref. |23|. In this approximation the partition function 
becomes 



exp 



knT 



I?(/iexp — 



1 



knT 



L/2 



dz 



L/2 



C fd(j> 
dz 



(2.7) 



where m — ai — Aa2- We shall consider the following ansatz for the form of 4>{z)] 



7z 



(z) = 0p(z) + — + (/)o. 



(2.8) 

where (t)p{z) is the component of (f){z) that satisfies periodic boundary conditions, i.e. (f>p{—L/2) = (f>p{L/2) and it 
has a spatial average of zero along the length of the molecules. 
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Having a period of L, (j)p{z) may be expressed as 

0p(z) = ^^&„exp — — . (2.9) 



The (t)o component of (t){z) is used for the calculation of rigid body fluctuations of the average relative azimuthal 
orientation between the two molecules and so does not depend on z. The component of 4>{z) proportional to 7 is 
included in order to account for free fluctuations of the ends by allowing 7 to vary, where 7 is the difference in the 
azimuthal angles between the ends of the molecules. We may recast Eq. (2.7) in the following way, using Eqs. (2.8) 
and (2.9), 



Z = exp (^--^^ j d<t)Q J djYi J ^ 




2 + 2 ^ 
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Lm 
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-1 L/2 

where G„ = (f (^)^ + m) and J„ = / dz(f)exp| 

The crucial step is to make the variable shift bn bn mjJnGn to decouple 5„ from 7. In the domain 
—L/2 < z < L/2, ^ may be written as a Fourier series and the effect of the variable shift is to adsorb this into 
4>p{z). At z = L/2, however, ^ cannot be expressed as a Fourier series. At this point there is a difference of 7 in the 
sum of the Fourier series representing ^ and the actual value of The torsional energy term, however, contains 
derivatives of which are sensitive to this difference. Therefore, when we make this shift we flnd that the integrand 
in our partition function still depends on 7. Thus this leads to a difference in the partition function for free boundary 
conditions as opposed to that assuming periodic boundary conditions [29l |. 



The partition function thus becomes 



Z = ZoZfZpexp[-^], (2.11) 



where Zq — dcjjQ exp ^— 2k^T ^'oj corresponds to the component of the partition function that depends only on the 
average azimuthal orientation of the molecules, Zf — J dj exp I — -^^^ coth (^) ) is the part of the partition function 



that takes into account the free rotations of the ends of the molecules , and Zp — Yl J dbn exp | — 51] '^'g " I 

is the component that takes into account all z-dependent azimuthal angle orientations in the case where the ends of 
the molecules share the same orientation. 

From Eq. (2.11) we may compute the free energy (Appendix B). Wc find that 

CL keT /A„\ keT ( ( 2L\\ 

F . -k.TO + ^0 + ^ + ^ m (^) + ^ m (1 - exp (2.12) 

Expanding the log of the last term when L is large, we retrieve Eq. (2.6). And so we find that this formulation is 
consistent with the exact one of Ref. 0|. We therefore are ready to tackle the many body problem of assemblies. 
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FIG. 1: Lattice labeling and relative positions of the lattice vectors u and v . 



III. ASSEMBLIES OF MOLECULES OF FINITE LENGTH: THE GAUSSIAN APPROXIMATION 

For an assembly, the energy can be written as a sum of the energies of each molecular pair over the entire lattice 
that defines the assembly as 



E\ 



C fd(b,i{z) 
dz 



^ (-I)"a„(i?2)cos {n{(i)ji{z) - 



+ (-I)"a„(i?i) (cos {n [c^ji{z) - (i>,^xi{z))) + cos {n {i>fl{z) - 0,7-i(z))))] . (3.1) 

Here, we have introduced two lattice vectors ui — jRiu and Vj = IRiv, where ii and v are unit vectors that describe 
the relative positions between molecules situated at the sites of the two dimensional lattice. These two vectors, as 
well as the convention used in labeling the lattice sites, are shown in Fig.l. 

In Eq. (3.1) we have included only the interactions between nearest neighbors and have assumed that all the DNA 
molecules are of the same length with the center of each molecule lying at z = js^l . Again we may truncate the 
series dX n — 2. In calculating the energy of the lattice, we have allowed for two separations, i?i and Ri, of the six 
nearest neighbors of a given molecule, since there may be distortions from the hexagonal lattice [iM l20l|. For this 
distorted or "rhombic" lattice, Ri corresponds to the distance to the neighboring molecule at an adjacent corner of 
the unit cell while R2 is the distance across the short diagonal of the rhombic cell. We may define a distortion angle lu 
(shown in Fig. 1) that characterizes the relationship between the two separations in Eq. (3.1), R2 = Rw/^, — 2cosci;. 
In the ground state, the amount of distortion is determined by minimization of with respect to the distortion 
angle uj. The relative strengths of the a coefficients determine the degree of this distortion, e.g. the oq term, which 
contributes a repulsive component to the overall force, has the effect of reducing the amount of rhombic distortion as 
on its own it favors a hexagonal state. For certain azimuthal configurations of the molecules at large densities, the 
distortion angle lo is greater than 60°. 

Before considering effects of thermal fluctuations, let us first describe the various ground states that the assembly 
has. When molecular separations are large, all the relative azimuthal orientations of the molecules are the same. 
This state may be referred to as the "ferromagnetic" state. For this state the lattice is hexagonal, Ri = i?2. Below 
a critical value of the separation of nearest neighbors, R\, the assembly adopts a new state where four of the six 
nearest neighbors about a specific molecule adopt a different azimuthal orientation to that molecule, as shown in Fig. 
2. This state then may be termed the "antiferromagnetic" state. This configuration of the azimuthal orientations of 
the molecules favors rhombic distortions, Ri 7^ i?2, described in the previous paragraph. The ground state azimuthal 
configuration for this configuration is characterized by the following equations 

{<j>ji{z) — (j>j^ii-i{z)) ~ (across the short diagonal of the rhombic unit cell) (3.2a) 

{(f)ji{z) — cf)j^ii[z)) = {(t)ji{z) — (j)ji^i{z)) = ip (adjacent corners of the unit cell) (3.2b) 

where cos(-0) = For even denser assemblies, a new critical value of average separation R^ is realized where the 

most energetically favorable state is the "Potts" state. In this state, the lattice returns to a hexagonal confirmation. 
As shown in Fig. 2, the azimuthal orientations of the molecules for this state may have one of three different values: 
(/)o, 00 + ^p, and (f>o + 2ipp where 



cos(V'p) = ifl+\/l + — ). (3.3) 



4 V V 0.2 
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FIG. 2: Schematic pictures of the azimuthal orientations of molecules in two of the ground state configurations. The first (a) 
corresponds to the antiferromagnetic state where the molecules adopt layers in which the molecules are in the same azimuthal 
alignment. The second (b) corresponds to the Potts state, where the azimuthal orienations of the molecules have one of three 
values as defined in the text. 



Upon introducing thermal fluctuations, we shall only consider the ferromagnetic and antiferromagnetic states. (For 
a discussion of the Potts state, which, again, only occurs at very large densities, at finite temperatures for rigid 
molecules see Ref. poj.*) First, for the antiferromagnetic state, the partition function in the Gaussian approximation 
may be written as 



Z — exp — 



knT 



n 



Vcj) {z) exp 



TT^AF 



(3.4) 



where 



L/2 



dz 



-L/2 



d(t)'M)\ TO2 



dz 



/.;,(z) - 0;._i,(z)) V (<^;,(z) - ^;,_i(z))' 

(3.5) 

and the ground state energy is E^^ = ao{R2) — ai(i?2) + a2(-R2) + 2ao(i?i) — 2ai(i?i) cos(V') + 2a2(i?i) cos(2?/'). The 
azimuthal orientations of the molecules at the lattice points are replaced with the following values: 



(j)jl{z) - (j)j-iiiz) = V + (t>'ji{z) - 0j_i;(2) 
(j)jl{z) ~ (j)ji-i{z) V + (t>'ji{z) ~ (t)'ji-iiz) 
(j)jl{z) ~ (/)j+y_l(z) = (l}'ji{z) - (jj'j+ii^iiz). 



(3.6) 



Here, 7712 = ai(i?2) — 4a2(i?2) and m2 = ai{Ri) cos{tp) — 4a2(i?i) cos(2-0). As in the previous section, we introduce 
the following ansatz, 



<i>',i{z)^r,^{z) + :^ 



(3.7) 



splitting up the azimuthal value along each molecule at each lattice site into a function (j)^i{z) that satisfies the periodic 
boundary conditions, (/)^;(— L/2) = (/)^j(L/2), linking the ends of the molecules and the correction to account for 
independent fluctuations of the two ends of each molecule. We may express (f'^iiz) in a similar form as Eq. (2.9) for 
the previous section (see Appendix C). The rigid body part of 4>ji{z) has been included in (j^iiz) and is the spatial 
average of 4>^i{z) along the molecule. The steps in calculating the free energy in the Gaussian approximation are 
similar to those given in the previous section (see Appendix C) with the end result being 
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ksT 
2(27r)2 

ksT 
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'C{x,y;ay/^L 



C{x,y,a) 
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ksTQas 



(3.8) 
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^here Al = ^ 2{a,iR,)cosW-4a,(R,)cos{2i,j ) ' ^2 = ^ 2[a,{R.,)-ia,^ ^ ^'(x, y; = COs(x)) + (1 - COs(y)) + a(l - 

cos(a: — y)) and a = (Ai/A2)^. For the ferromagnetic state we simply set '0 = and Ri = i?2- 

Additionally, we may compute the correlation functions associated with the thermally induced fluctuations of the 
azimutathal orientation of the molecules. Calculation of these are particularly useful since they are reflected in the 
intensity profiles of x-ray diffraction patterns of columnar assemblies [lTI | . This azimuthal correlation function within 
the Gaussian approximation is defined as 

(exp {^n (0;,(z) - 4,(z'))))o = ^xp (-^) 1 11 / ^ " exp (-^^f^) • (3-9) 

This may be written as (summarized in Appendix D) 

(exp(^n(0;,(z)-0;.,,(z'))))o - e^v{n' 

= e^p{n^GU-fJ-l',z,z')). (3.10) 
G{j — j' , I — I' , z, z') may be split up into terms corresponding to three contributions 

G{s,r,z,z') = Go(s,r) + Gp{s,r,z- z') + GAp{s,r, z, z'), (3-11) 
where, s = j — j' and r = I — I' . 

The first contribution takes into account rigid body fluctuations of the azimuthal orientation of the DNA molecules. 
It has the following form for large separations, i.e. s, r ^ 1 (for a general expression see Appendix D) 

where = Go(l,0) = Go(0, 1) = 2|^^arcsin ^-^=1=^ The flrst term in Eq. (3.12) grows with increasing 

separation because of there is no long range order for such a 2-D system. The second contribution of Eq. (3.9) 
represents the torsional (z-dependent) azimuthal fluctuations assuming periodic boundary conditions. Again, when 
the separations are large, 

Gp{j^y,i-i',z-z')^-^l 



2C 

— TT — 7r 



The last contribution is the correction due to allowing both ends of each molecule to fluctuate independently. This 
takes the form for large separations of 

""'^^''^ ' ^ = ^ (^ I'^l '' Gi.,y;ay/^ [ A^ j [' " [ A^ j j 

— TT —IT ^ / \ \ / / 

exp nA£i^E^\ „ exp (^-^^:igi^));3.14) 



z'\C{x,y;ay/^\ ( ~ \z'\G{x,y-af/^ 



A^ j-^'^PV 

And so the correlation function at large separations can be written as 

(exp [^n (0;,(z) - 4.(z'))))o - - j' f + H - I'f + ^^^^^^(Tt|^) ^ 

X exp {-n^ (Al + )) exp {~n^G^p{z, z')) , (3.15) 

where 7 = n^ksTXi/ {TrLG\/l + 2a). As the separations increase, the azimuthal correlation is reduced as one might 
expect. This, again, reflects the fact that in such a 2-D system there is no long-range order. 
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IV. ASSEMBLIES OF MOLECULES OF FINITE LENGTH: THE SELF-CONSISTENT 

APPROXIMATION. 



Near the point of frustration, i.e. the location of the transition between the distorted "antiferromagnetic" state 
and the hexagonal "ferromagnetic" state, the fluctuations become quite large so that the Gaussian approximation is 
no longer valid. Thus we must extend the calculation by using a self-consistent approximation, namely, a Hartree 
approximation. Unfortunately, we were unable to find a Hartree approximation that takes into account freely 
fluctuating ends of the molecules, but it is possible to treat this contribution as a small correction to the Hartree 
approximation for periodic boundary conditions. In such an approximation it is possible to determine how torsionally 
softening or changing a molecule's length alters the value of Rl{T), the separation at which the phase transition 
occurs. Furthermore, the correlation functions of the previous section also can be calculated within this approximation. 

The self-consistent approximation first entails carrying out a series expansion of Eq. (3.1) in powers of (j)' where 
anharmonic terms, terms of 0{4>'^) or greater, are treated as perturbations to the Gaussian approximation. The value 
of tp is determined through the requirement that 



(4.1) 



(Details of such a perturbation expansion are given in Appendix E) Terms for the correlation function and the free 
energy, where periodic boundary conditions are assumed, may be re-summed in a similar fashion to the Hartree 
approximation calculations of Refs. ,20t. .23i | . Following similar steps (given in Appendix F), we obtain the following 
Hartree result for the free energy per molecule 
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where au = I ^2)^ ^ ^'^'^ ^1 ^'^'^ ^2 satisfy the following transcendental equations 



Af = 
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il)H satisfies the following relation 

COS (1/)^ ) = — exp 



ai /3Af 



4a2 



V4Ap^' V2Af 



L 



(4.2) 



(4.3) 



(4.4) 



The functions xi (^T^^^'^h^ X2 (^J^^'^h^ given in Appendix E. The free energy is then minimized with 
respect to ui to determine the degree of distortion of the hexagonal lattice for the antiferromagnetic state. The 
ferromagnetic state is simply obtained by setting -0 = so that Ri — R2 and, thus, Af = Af = \h and a = 1. 



It is quite interesting to look at the Hartree result for the ferromagnetic state for infinitely long molecules, which 
is very similar to the result given in Ref. '23l for a pair of molecules, since it provides some insight into the nature 
of the many body effects of the assembly. There are two differences: first of all, Ap is replaced by ApC~J, where 
Cmb is a constant resulting from the many body effects of the assembly; and secondly, an overall factor of three that 
multiplies the free energy is included to account for the coordination number of nearest neighbors about a molecule. 
We calculate that Cmb — 0.560 so we find many body effects suppress thermal fluctuations. This is not surprising 
considering that the many neighbors in the assembly should increase the effective interaction between molecules. 
Therefore, thermal fluctuations of the same magnitude in the assembly cost more energy than that for a single pair 
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FIG. 3: Location of the antiferromagnetic-ferromagnetic transition for molecules of different lengths and torsional flexibilities. 
At large densities (small interaxial spacings), the antiferromagnetic state is the lowest energy state whereas the ferromagnetic 
state is the favored state for more dilute assemblies. Diamonds (o) correspond to perfectly rigid molecules (C = oo), stars 
(★) to molecules with a torsional modulus of C = 3.0 x 10~^^ erg-cm, circles (o) to molecules with C — 1.0 x lO"'^^ erg-cm, 
triangles (A) to molecules with C — 3.0 x 10~^* erg-cm, and squares (□) to molecules with C = 3.0 x 10~^^ erg-cm. Results 
are shown for molecules with 70% charge compensation and a 70/30 major/minor groove charge distribution in a solution with 
an inverse Debye screening length of about 7A. 



of molecules. 

We may also show that the results above conform to that of the formulation of rigid molecules |23|. Changing 
variables to and where Af'^ — \/ CLj (2J^) and — \/ CLj (2 J|^), we may rewrite Eq. (4.3) in terms 
of and ■ Upon takin g th e limit L ^ 0, it is fairly straightforward to recover the Hartree results for the 
configurational states of Ref. |23|. Taking the limit of these equations, i .e. the molecules are treated as completely 
rigid, also yields this same result. This is fully consistent with the physics we might expect. As we shorten the length 
of the molecule, the energy cost of torsional fluctuations begins to counter the reduction in free energy due to the 
entropy gain associated with them. And so with decreasing length, these fluctuations become suppressed and rigid 
body fluctuations start to dominate. 

We first examine how the phase boundary between the ferromagnetic hexagonal state and antiferromagnetic rhomic 
state changes with the length and flexibilty of the molecules assuming periodic boundary conditions. Equating the 
free energy of the antiferromagnetic state with that of the ferromagnetic state, Eqs. (4.2)-(4.4), we can determine the 
location of the transition between the two states. The location of this transition in terms of the average interaxial 
spacing between molecules in the assembly as a function of molecular length is shown in Fig. 3 for molecules 
with varying torsional flexibilities. Increasing the flexibility of the molecule results in the transition moving to 
larger densities (smaller interaxial spacings) as one might expect. With increasing flexibility, there is an increase 
in thermally induced torsional fluctuations and so the location of the transition as the flexibility is increased moves 
in the same direction as when the molecule is shortened, where thermally induced rigid-body fluctuations become 
larger. Also shown, in Fig. 4, is the level of distortion (value of the rhombic angle ) from the hexagonal lattice as 
a function of the average interaxial spacing for the antiferromagnetic state. As seen in the plot, greater flexibility 
reduces the amount of distortion. 

Modifications of the results of the self-consistent Hartree approximation when incorporating corrections for 
independently rotating ends of the molecules are quite involved and are left to the appendices. When including these 
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FIG. 4: The level of distortion (uj) from the hexagonal lattice for the antiferromagnetic state as a function of the average 
interaxial separation of the molecules in the columnar assemblies. The distortion is shown for the ground state configuration 
(no spin fiuctuations) , for a rigid molecule of length 500A with spin fluctuations, and for a flexible molecule (C = 3.0 x 10~^^ 
erg-cm) with torsional fluctuations of the same length. 

corrections, however, we find that the locations of the transition only shift by hundredths of angstroms for all the 
cases shown in Fig. 3. And so the formulation assuming periodic boundary conditions, as long as the molecules are 
long enough, adequately describes the system. 

We may also calculate correlation functions within the Hartree approximation for periodic boundary conditions. 
Here, we may use the results of the previous section, but we now replace Ai with and a with an = {^1/^2)^- 
We then may compute and within this approximation. A plot of these quantities is shown in Fig. 5 for a 
molecule with a realistic value for its torsional modulus. The correction associated with fluctuations of the free ends 
of the molecules only results in small changes in the value of A^. (The formulation of which is again left to the 
appendices.) On this plot, we also show the variation of the relative azimuthal angle due to rigid body fluctuations, 
A^ , for a molecule of the same length. These results are consistent with fits of x-ray diffraction data of hydrated 
DNA assemblies [ll| . 

Lastly, we show in Appendix H the high temperature expansion for the azimuthally disordered or Kosterlitz- 
Thouless vortex phase, mentioned in the introduction, that includes torsional fluctuations. Granted, in the case of 
real DNA it is very unlikely that this situation is reaUzed M- At larger separations where thermal fluctuations 
dominate intermolecular azimuthal interactions, the DNA assembly is likely to lose columnar ordering before such a 
transition in the azimuthal configuration of the columnar assembly occurs. Because of this reservation, this work has 
not been included in the main text. However, in other helical bio-molecule assemblies where azimuthal interactions 
might be weaker but spatial interactions are strong, a Kosterlitz-Thouless like transition could possibly occur within 
a smectic layer with short range two-dimensional ordering. 
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FIG. 5: The contributions of torsional fluctuations (fiUed sqares) and rigid body fluctuations (open circles) to the 
correlation function of the relative azimuthal angles between molecules in an assembly as a function of the average molecular 
separation. Again the break at 25. lA corresponds to the transition from the antiferromagnetic state (at smaller molecular 
separations) to the ferromagnetic state. Values are shown for an assembly of molecules that are 500A long and have torsional 
moduU of 3.0 x 10~^® erg-cm. 



V. DISCUSSION AND CONCLUSIONS 



We have developed a method for calculating the contribution to the free energy due to azimuthal fluctuations 
for torsionally flexible molecules of finite length. Using this method we can determine the influence of azimuthal 
flexibility of molecules within a DNA columnar assembly on the position of the phase transition between the 
antiferromagnetic state where the hexagonal lattice is distorted and the ferromagnetic state for various molecular 
lengths. For molecules with torsional moduli of C = 3 x 10~^^erg — cm, the approximate value deduced from 
experiments [3lj |. torsional flexibility reduces the separation at which this occurs compared to that found for the 
case of rigid molecules. We see that the effects of torsional fluctuations are more pronounced for assemblies of long 
molecules, whereas the results for short molecules are little different from those rigid body calculations of Ref. |20|. 
thereby showing that shortening the molecule is nearly equivalent to making the molecule less flexible. 

We found that torsional fluctuations considerably reduce the degree of distortion, for long molecules in the 
antiferromagnetic state, as compared to that found for calculations for rigid molecules of the same length. Since 
flexibility allows for readjustment of the azimuthal coordinate along the length of the molecules in the lattice, less 
distortion is required to minimize the free energy of the assembly due to the antiferromagnetic coup ling of the 
pair potential. Hence, as compared to ground state ^3 or finite-temperature rigid body calculations |2Clj |. X-ray 
diffraction patterns must have even better resolution to pick out this distortion 0|. 

We were also able to calculate the asymptotic form of the correlation function {exp {in {(j)ji{z) — (pj'i' {z')))) 
including both rigid body and torsional fluctuations. This term is found in the formulation of the intensity of x-ray 
diffraction patterns of assemblies |llj | and is therefore a direct means by which one can determine the extent of these 
fluctuations. Our results for the parameters chosen agree quite well with the fits made to x-ray diffraction data of 
hydrated DNA assemblies. Conversely, comparing these fits to the correlation function formulation may reveal the 
strength of the azimuthal interaction, which in turn may say something about the charge distribution on the DNA 
molecules in the assembly. 

Although we have discovered a number of effects associated with allowing for torsional flexibility in columnar DNA 
assemblies, this analysis, however, is limited by the form of the DNA-DNA interaction we have employed. The 
effects of nonlocal polarizability of the water in the narrow interstitial regions between the DNA could alter the results 
[lEHHIlS' Also, the potential assumes a given charge distribution of readsorbed counterions that changes very little 
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with density. Once these effects have been incorporated into the pair potential the same analysis developed here 
may be used. But even upon considering these effects, the qualitative aspects of this analysis should remain unchanged. 

One may also need to consider the effect of sequence dependent distortions from the ideal double helical structure 
seen in real DNA ,3^ ,3^ ,3^ 35] in the analysis of the assemblies. This may well be included in calculations for 
assemblies of non- homologous DNA by employing ensemble averaging J^- Such averaging is easiest in an assembly 
made up of with several different DNA sequences. Again, as in Ref. l23'|, these effects may be incorporated into a 
combined persistence length that incorporates both distortions arising from thermal fluctuations and from different 
base pair sequences of neighboring molecules. For assemblies containing DNA sequences with only a few different 
texts (and thus different distortions from an ideal helix), the interaction is much more difficult to treat and has yet 
been considered. Nevertheless, in the case of homologous DNA assemblies, distortions from an ideal helix make lit- 
tle difference as the electrostatic interaction depends only on the relative azimuthal orientations between the molecules. 

Finally, we should point out that this is by no means a complete statistical mechanical picture of a columnar 
assembly. For example, we have not included lattice vibrations, fluctuations in the ^-position of the center of each 
molecule, bending fluctuations, nor fluctuations due to other geometrical distortions of each molecule. However, 
fluctuations in the azimuthal degrees of freedom are likely the most important feature for determining the structure of 
the mesophases of columnar assemblies as seen in the fits of x-ray diffraction patterns |ll|. Furthermore, coupling of 
the azimuthal degrees of freedom to other geometrical distortions, which are likely to be small, may be incorporated 
into our model. 
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Appendix A. Details of the "Quantum Mechanical" calculation of the free energy at finite length. 



We start with equation. (2.5) of the text, namely 



F = -kBTln J2 I cxp ( -— ) {E I 0-) - kBTlnQ. 



(Al) 



We will now, using Eq. (Dl), derive expressions for finite size correction terms for the Gaussian approximation. 
But first we will need to derive some general results. 

Now each Eigen state of Eq. (2.3) is characterised by two numbers; the band number 6 > 0, an integer, and the 
wave number in the band 1/2 < k < 1/2, a continuous variable. Each Eigen state should then be written as 



V'fo,/c(0) = (0 I Eb,k) = fim 1— G„,b,k exp {i{n + k)(t)), 



(A2) 



where the Gn.b,k satisfy 
'{n + kf 



El 



b,k 



Ap fl2Ap 
Gn,b,k — -^^YJ [Gn-l.b,k + Gn+l.b,k\ + -Z ^2 [Gn-2,b,k + Gn+2.b,k\ — 0. 



2\l 



(A3) 



Here, Np is the number of periods the potential has in (/)-space and is proportional to the system size: the range of 
values which cf) is allowed to take. As cf) will fiuctuate between —oo and oo, except at 0+, we take Np to be infinite. 
It is easy, then, to show that the following is true 



27r 
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ipb,k{(l))d(j) = lim -= 
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^k,oGob,0- 



(A4) 
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Using Eq. (D4) we rewrite F as 

F^-fcsTln^lGo.Ml'expj^-:^^ -fesTlnB. (A5) 

In the Gaussian approximation we assume that 4'b,o can be written as a hnear supposition of simple harmonic oscillator 
Eigen states 

V'M- J2 i^b(^^^(<^-27^,))expf^^^^^^-^Y (A6) 

where the Hb{x) are the Hermitc polynomials and the Cb are the normalization constants for each state. In the 
Gaussian approximation 



= I "r| I cos 00 — ( — -Tj- 1 cos 2(f}Q and 



2 



Eb = [- + b]u;- {^] cos(</)o) + ^ cos(20o). (A7) 



Assuming negligible overlap between the wave- functions in the superposition; we find on retaining the first two terms 
in the sum in Eq. (D5) and expanding out the logarithm and writing lo — Xp/X we find 

^ ^ CL keT / (E2-Eo)L\ ksT , f X\ , , ^ , 

We should point out that this result differs from (D9) of as the previous version a couple of mistakes that have 
been corrected. Substituting for Eq and E2 we so obtain Eq. (2.6) of the text. 

Appendix B. Calculation of the free energy for the DNA pair by Gaussian integration. 

We may readily perform the Gaussian integrations in Eq.(2.10) of the text 

Z = Z„ZfZpexp(~^Y (Bl) 



where now 

Zo 



Zp^\{ {2^kBTG^f'\ (B2) 



Using the definition of free energy F = — fc^TlnZ we are then able to write the free energy as 

where we have adsorbed terms that dont depend on A and are not important to our analysis into O. The first term 
may be evaluated by means of a trick 

d ' - ■ ^ ^ .... s -1 



dm 
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We then by integrate up (B3) and so arrive at 

F = ..Tin (si„h (£) ) + M (I c„.h (i) ) + tere, (B5) 

where extra constants, such as constants of integration have been adsorbed into O. On rearrangement we arrive at 
(2.11) in the text. 



Appendix C. Calculation of the Free Energy for Gaussian fluctuations in the eissembly. 



We start our analysis by substituting (3.7) into (3.5) of the text and so obtain the following expression for [(/>'] 



dz 



+ 



2L2 
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We then introduce the following lattice Fourier transforms 



(CI) 



(C2) 



fe„fe„ 



where A is total area of the lattice and ro is the lattice spacing. We have chosen reciprocal lattice vectors corresponding 
to the rhombic Bravais lattice defined by u and v see Fig. 1 of the text. Then fc„ and ky take on values which lie 
within the first Brillouin zone for a rhombic lattice (e.g. — 7r/ro < ku < 7r/ro). We are then able to write in the limit 
00 



-TT — TT 
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(F{x,y,z) + 
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2L2 
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where x = r^ku and y = r^ky. We then may express <?^(a;, y, z) through Fourier Series 

/ 2'Kinz\ 
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and so re-express (C3) 
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C Lm^, 
■j^ + —C{x,y;a 



l{x,y)j{-x,-y) 
2L2 



where Gn{x, y) = (c {^f + 2mC{x, y; a)) \ J„ = (1 - Snfi), C{x, y; a) = {1 - cos{x)) + (1 - cos(y)) + 

a{l — cos{x — y)) and a = m2/mi. On making the variable shift bn{x, y) — > hn{x, y) — mi'y{x, y)C{x, y; a)JnGn{x, y) 
we may write this as 
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We find that on further manipulation we may write 



f]AF \ 

Z = ZfZpexp[-^\, (C7) 



where 
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and S{x,y) = x" c(x y-ay/^ tanh ^^C'(a;,t/;a)^/^y We may evaluate the path integrals giving 

lim TT (27rs(-^,-!^-^ 



1/2 



The product iV = N^Ny is the number of molecules in the assembly which we have taken to be infinite as ^ — > oo. 
From (C8) we find the following result for the free energy per molecule 
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+A;BTe„,6, (C9) 

where unimportant terms that do not depend on mi or a have been adsorbed into Qasb- Using the trick illustrated 
by Eq. (B3) It is possible to write 

F=^. Ux [dyln(^ (cix,y;a)Y\^,Jcix,y;a)^f^L 



2(27r)2 J i ^ I A \ C{x,y;l) / I 2A 

TT TT / ^ 

+ ^ Jdxj dy\n \^^J-—smh(^-^C{x,y,a)U+kBTe,sb, (CIO) 



-TT — TT 



where additional constants have been adsorbed into ©os6- 

Appendix D. Correlation functions for Gaussian fluctuations in the assembly. 

It is first useful to consider the following correlation function 

Gij-f,i- i',z,z') = (<^;, - 1 {M^f)o - \ (<^/<'(^')')o ' (Dl) 

where in general for any quantity A\(^ 

= exp (-^) In/ T^'^^MAm exp (-^^^) , (D2) 
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using the ansatz (3.7) we may write (Dl) as 
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This may be rewritten as 
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We then make the variable shift 6„(a;, y) bn{x, y) — mi'y{x, y)C{x, y; a)JnGn{x, y) and obtain the following 
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In these expressions we are able to perform the sum and we also separate out the rigid body contribution by writ- 
ing Gp{j,l,z) = Gp{j,l,z) + Go{j,l) where Go{j,l) is the part of G{j,l, z, z') that arises purely from rigid body 
fluctuations. We obtain the following results 
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Now from these expressions let us find their behavior when l,j ^ 1. For both GAp{j,l, z, z') and Gp{j,l,z) this 
is simple. As j and / increase, the terms in the integrands of GapU, I, z, z') and Gp{j, I, z) that depend on j and I 
oscillate more rapidly, so that their contribution gets smaller and smaller. So when l,j^ 1 we may neglect terms 
that depend on j and I in both integrands, so arriving at (3.13) and (3.14) of the text. The term Go{j,l) is trickier 
as neglecting terms that depend on j and / leads to a logarithmic divergence. It is, however, possible to do one of the 
integrations in (D8) (by contour integration around a unit circle) which leads to the following result (for j > 0) 
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Now when the integral is dominated by small values of y. We may expand out the integrand for small y, and 

so obtain 
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where r = 1 + aj/{a + 1). 



Here, it is convenient to subtract off Go(l,0) as we shall see. We may evaluate (Dll) by differentiating (Dll), 
which enables us to perform the resulting integrals, and so we obtain the following differential equations, 



2kBT 



(1+a) J 



\djj^ {2n)Lmi{l + a)j^±^f + r^' 
dl _ 2kBT r 

~ {2Tr)LmiVT+2^^g^j^ + ' 



(D12) 



with the boundary condition that /(1,0) = 0. By solving (D12), subject to the boundary condition, substituting for 
r, we find the following asymptotic form 



knT 



(27r)LmiVl + 2a 



In f + r + 



2ajl 



(D13) 
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and so we are able to obtain (3.12) of the text. 

We find that able to show that on applying (3.5), (D2) and (C6) and making the shift bn{x,y) — > bn{x,y) — 
mi'y{x, y)C{x, y; a)JnGn{x, y) that we may write 

(exp {m (</.;,(z) - <^'M^')))), = |exp (-^) \{ J Vhn{x,y) J Vj{x,y) 



e^pi^jdxj dy ^^^^y (b„{x,y) +j{x,y)Jn{x, 

y)) (D14) 

\ ^ — TT — TT 

expf--i- Jdx 7(x,y)7(-^ -2/) _ 1 ^ J bUx,y)b.^i-x,-y) 

^1 {27Tf J J 2S{x,y) (27r)2^; J 



-TT — TT — TT — TT 



2Gn{x,y)kBT 

where Jn{x,y) = J„(l — 2mC{x,y;a)Gn{x,y)). Now on completing the square and using (C7) we may show (3.9). 



Appendix E. Perturbation Theory. 

The interaction energy may be expanded out to beyond quadratic order in (j)' according the prescription given by 
equations (3.6) and (4.2). 



Z = exp 



n /-^'(^)exP (-^i^MMM) . (El) 



El 



The first two terms, Eq^ and E-^^[(f)] are given in the text. The next term, E-^^[(p] comes from linear terms, which 
no-longer vanish as (3.2) is no-longer satisfied. This term takes the form 

L/2 

Et^W]=Y. J (ai sin (V) -2a2 sin (2V>))((<^;,(^) (^)) + (<^;7(^)-0;,-i (-2))). (E2) 

•^^ -L/2 

The last term may be split into two pieces E-^f^ [(j)'] = E^jj [(j)'] + E^jj [(j)'] where we have 

L/2 

E^2lW]=Lj2f2 j ^^(^x{(«i(-ir"'+«2(-4)")(</.;.,,(^)-0;.+i,,_i(^))'"+ (E3) 



jl n=2_ 



(2n)! 

-L/2 



+ (ai(-l)"-i cos(V') + cos(2V')a2(-4)") (((<^;,,(^) - 0;_i,;(z))2" + (<^;. ,(^) - cl>'.^i_,{z)f-))} 



L/2 

E^ahW]=lJ2Y1 / (2n - 1)! ^^"^("1)"'' + "2 sin2V>(-4)"/2) 

((<^;.,(^) - <^;_i,;(^))^"-^ + (0;.,K^) - . 

Here, and throughout the appendices, we adopt the convention a„ = a„(i?i) and a„ = a„(i?2)- For the moment 
let us truncate these scries at n = 2 terms. We then use the ansatz (3.7) and perform the lattice Fourier transforms 
given by (C2) as well (C4).We then may write 

E^^[(j)'] = (ai sin (ip) - 2a2 sin (2^')) Vl lim [(2 - exp(ia;) - exp(iy)) bo{x, y)], (E4) 
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TT TT TT 7T TT IT 



-IT — TT TT — TT — TT TT 



+92F^'^\x,x',x",y,y',y")^ {hn{x,y)bn'{x' ,y') bn"{x" ,y")b-n-n'-n"{-x - x' - x",-y -y' - y") 

+46„(a;, y)bn'{x', y')bn"{x", y")^{-x - x' - x" , -y-y' - y")J-n-n'-n"{-x - x' - x", -y-y' - y") 

+Qbn{x, y)bn' {x', y'hix", y"h{-x - x' - x", -y-y'- y")Jn" {x", y") 
J-n-n'-n"{-x - x' - x" , -y - y' - y") + 46„(a;, y)^{x' , y')'y{x", y")^{-x - x' - x" , -y-y' - y") 
Jn- {x', y')Jn" {x", y")J-n-n'-n" {-x - x - x" , -y - y' - y") + 'y{x, y)-f{x' , y')j{x", y") 
■y{-x -x' - x", -y-y' - y")Jn{x,y)Jn'{x' , y')Jn"{x", y")J-n-n'-n"{-x -x' - x", -y-y' - y"), 



TT TT TT TT 



4'irM = -^^;^^E Jd^J d^'J J dy'9^'^F^'\x,x',y,y') 



TT — TT — TT — TT 



{bn{x, y)bn' (x', y') 6_„_„' {-x - x', -y - y') + 37(0;, y)Jn{x, y)bn' {x, y)b-n-„' {-x -x' ,-y- y') 
+87(0;, t/)7(a;', y')Jn{x, y)Jn' {x' , y')b-n-n' {-x - x' , -y - y') + 7(2;, y)-i{x' , y'h{~x ~x',-y- y') 

Jn{x,y)Jn'{x' ,y')J-n-n'{-x - x' ,-y - y'), 



where 



F^^\x,x',y,y') = (1 - cxp(-v;.t))(1 - exp(-ix'))(l - exp(i(x + x'))) (E5) 
+ (1 - exp(-it/'))(l - exp(-iy'))(l - exp(i(y + y'))), 

f[^\x,x' ,x' ,y,y' ,y') = (1 - exp(-ia;))(l - cxp(-ia:;'))(l - exp(-ia;"))(l - exp(i(x + x' + x"))) 
+ (1 - exp(-i2/))(l - exp(-i2/'))(l " exp(-zy"))(l - exp(i(y + y' + y"))), 

F^^\x, x', x', y, y', y') = (1 - exp(-z(x - y)))(l - exp(-z(x' - y')))(l " exp(-z(x" - y"))) 

(1 - exp(i(a; -y + x' -y' + x" - y"))), 

g^^^ = aisin^/; — 802 sin 2^, g^^ = aicosf/' — 16a2 cos 2^/1 and g^^ = di — 16d2. Let us calculate the leading order 
correction to the Gaussian result for 

^nix,y) = exp (^-^^^ ^YiJ 'Dbn{x,y) J Vj{x,y)bn{x,y)b-n{-x, -y) 



when we neglect E^^[(j)] and set = 0, in the Gaussian approximation, we find that G^{x,y) = kBTGn{x,y). 

We derive the leading order correction by expanding out E^jj [0] 



Gni^^y) = kBTGn{x,y) + ex.p [ ^Yl J 'Dbn{x,y) / 2?7(x, y)6„(x, y)&_„(-x, -j/) 



F^ahW\ Ef^[(j)\ do not contribute directly to the correction, but do contribute to tp. Using standard procedures 
we may do the path integrations and so obtain the following result 

G^{x, y) = keTGnix, y) + {kBTfGn{x, yf{\{x, y) + Aj(x, y)), (E8) 



21 




FIG. 6: Diagrammatic representations that contain the leading order corrections to the Gaussian results for the full correlation 
functions Gn{x,y) and S^{x,y) defined in the text. For G^{x,y) the first graph from the left represents the Gaussian result, 
the middle graph the correction from Ap , and the last graph the correction from A/. Similarly, for S^{x,y) the first graph 
from the left represents the Gaussian result, the middle graph the correction from Ep and the last graph the correction from 



where 



TT 7T 



^^=2(2^^ / "^^'Z dy'Gn{x',y')(9[^^Fi^\x,-x,x',y,-y,y')+gi^^F^''\x,-x,x',y,-y,y')), (E9) 



-TT —IT 



and 



TT TT 



A/ = 2(^2nrL{kBT) ^ J '^'^'Z '^y'^^'^'^y'^ {g['^F^'\x,-x,x',y,-y,y')+gi^^F^^\x,-x,x',y,-y,y')) 



-TT — TT 



J-n'{-x',-y')Jn'{x',y'). (ElO) 

The first term represents the correction that comes purely from periodic boundary conditions and second term 
is the contribution from allowing both ends of the molecule to fluctuate independently of each other. Both these 
corrections may be represented diagrammatically in Fig. 6. 



We may perform the sums which then gives us 



A, 



2C 



sf'xi (^'"^^) [(2-2cos(a;)) + (2-2cos(2/))]+5^")x2 (^'"^^) (2 - 2 cos(a: - y) 



(Ell) 



L/2 Ljl 

A/ = j dx' I dy' {g['^ [(2 - 2 cos(a;)) + (2 - 2 cos{y))] (1 - cos(a;')) 



-L/2 -L/2 



+gi^'> [(2 - 2 cos(x - y)] (1 - cos(x' - y'))] } 



Ai 



2L C{x',y';a) 



fsinh ( L^i^l^yli^] cosh ( ^^^<y'-^-y''^ 

4C{x',y';ay/^ \^ \^ 2Ai ^ \^ 2Ai 



The functions xi (^2^r,o:H^ and X2 (^2^r,o:H^ are given by 

TT TT 

( = Z Z iw'-^—' J Cix, y; a) V2 ' 



-TT — TT 

TT TT 



-TT — TT 



22 



We may also calculate the leading order correction to the Gaussian result for 

ksT 



exp - 



We find that 



S^{x, y) = S{x, y) + S{x, y)'(Sp(a;, y) + Ef{x, y)), 



where 



" 2(2^)2^ / J dy'Gn{x',y') (^g[*^Fi{x,-x,x',y,-y,y') + g^2*^F2{x,-x,x',y,-y,y'jj 

Jn'{x,y)J-n'{x,y), 



-TT — TT 



(E13) 



(E14) 



(E15) 



TT TT 



(fcBT)L ^ / '^^'Z '^V'^^^'^y') [gf' Fi{x,-x,x' ,y,-y,y') + gf' F2{x,-x,x' ,y,-y,y')^ 



-TT — TT 



{x' ,y')Jn'{x' ,y' {x,y)- 
In Fig. 6, we show how the corrections are represented diagrammatically. On evaluation of all sums we find 



A? 



TT TT 



2(277)2(7 / / "^y' (^i''^^^ - 2cos(x)) + (2 - 2cos(y)))(l - cos(x')) 

— TT — TT 

+52'^H2 - 2cos(a; - y))(l - cos(a;' - y'))) -f^ (x,y,x',y'; 



(E16) 



-2' 



(27r)^ 



where 



7T TT 

T J dx' J dy' (gf\{2 -2cos{x)) + (2 - 2cos(2/)))(l - cos(a;')) + 5^^^(2 - 2cos(a; - y)){l - cos{x' - y'))) 

— IT — TT 



TT f / / Al 

-ff [x,y,x ,y; — , 



1 ^^^jLC(.W; 

C(x',2/';a)i/2 \^ 2Ai ^ 



a 



1 1 

" (7(a:',y';a)i/2 [4 



1 

_C'(x',y';a) - ^(x,^; 



. 2Ai 



y'; a) V2 coth ^ _ ^.^^^ ^. „)i/2 ^^^^^ |^ 

H I """^ 



a) 



,(lC{x\ 

r°«M ^ 

-—1^ coth { f sinh f ^ 



Ai 

" iC(a;,y;a)V2 



E17) 
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r'' = X + ^ + 



9 



T 



F = + 




FIG. 7: Diagrammatic representations that contain the leading order corrections to the Gaussian results for Tf{x,y) and the 
free energy F defined in the text. For Ff (a;,jy) the first graph from the left represents the Gaussian result, the middle graph 
the correction from Ap, and the last graph the correction from A/. For F the first graph from the left represents the Gaussian 
result.. The next graph along is Gi, that arises purely from periodic boundary conditions. Third from the left is G2 and the 
last is G3. These terms are corrections due to free boundary conditions. 

Now let us determine tp. We do this through the requirement that 

lim {(1 - exp(-iy)) {bn{x,y))} + lim {(1 - exp(-%))J„ {j{x,y))} | . (E18) 

Again, in the Gaussian approximation we may neglect E^H[(t>]- Equation (E18) is then satisfied by requiring that 

= ri = -lim {(2-2cosy)Go(a;,t/)}(aisinV'o -2a2sin2V'o), (E19) 

X— >0 

which implies that '^0 = or ^0 = arccos (ai/(4a2)), where tjjo is V for the Gaussian approximation, for Fi ti vanish. 
When we go beyond the Gaussian approximation to leading order in perturbation theory we find that Fi becomes 

Ff = — lim 1(2 — 2cosy)Go{x,y)} (oi sin^' — 2a2sm2ip) 

a:— >0 

-lim {(2- 2cos2/)Go(a;,y)} (aisinV'- 8a2sin2V') (ri +r'.) , (E20) 

X— >0 ^ ■' 

y^O 

where 

7r 7r 

r; = ^|f^E / J dy'{2-2cosy')Gn{x',y'), 

TT TT 

= ^j'^'^'j dy'S{x',y'){2 -2cos{y'))Ux',y')J-n{-x',-y'). (E21) 



-TT —TV 



Again, the first term, F^, arises purely from periodic boundary conditions and Fj is the contribution from allow- 
ing the ends to fluctuate independently of each other. Both these terms may be represented diagrammatically in Fig.7. 

We may compute the leading order correction to ip setting Fi = and by writing ijj = i/jo + i/j' and expanding out 
for small ijj' . We may also perform the summations in both F^ and Fj. We find that 



, / _ (aisinV>o-8a2 sin2V>o) , . 

^ - "(aicosVo-4a2Cos2V^o) ^ ^ + ' ^ ' 
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where now 



r' = ^1 

•'^ ~ (27r)24ApL 



7r TT 



dx'^ - cos(y')) 



-TT — TT 



Cix',y';a) 2Ai C(a;', y'; a)i/2 



(E23) 



2Ai 



2Ai 



We may also look at corrections to the free energy. We may write F = Fq — ksTAF where Fq is F evaluated in 
the Gaussian approximation (c.f (CIO)) and AF is the correction. We find for AF 



Gi = 



1 k^T 



AF = Gi + G2 + G3, 



TT TT TT TT 



(E24) 



j J J dx' J dy'Gn{x,y)Gn'{x',y') 



-TT — TT — TT — TT 



8 L(27r)4 

(gi^^F^"^^ (x, -X, x', y, -y, y'jj + gi^^F^'^\x, -x, x' , y, -y, y')) 



r - ^ ^ 

^ 4 L(27r)' 



TT TT TT TT 



^ I /^^ y"'^a;'y dy'S'(a;,y)G„'(x',y')'^(2;,y)-^-n(a;,y) 



-TT — TT — TT — TT 



{9i^^Fi^\x,-x,x',y,-y,y') + g^^^ F^'^\x,-x,x' ,y,-y,y')j , 



71 TV n TV 



^3 = ^ fcBTL(27r)4 ^ J J J J '^y''^(^'2/)'^(^''2/')-^n(a;,2/)4'(a;,i/)J„"(a;',y') 

J_„_„'_„" (x', y') (si'^^-P'i**^ (x, -X, x', y, -y, y') + g'i^ F^^^ (x, -x, x', y, -y, j/') 

These corrections are represented graphically in Fig. 7. The sums may, again, be evaluated yielding the following 

expressions 



L ksT ( Ai 



8 (27r)4 V C 



2g[\, ( 



(4). 



L 



(E25) 



^ L ksT fXi 



TT TT TT TT 



2 (27r)4 V^) I I I I 

— TT — TT — TT — TT 

- ™(x - - - „'))) _-L_ „„.h 

^^-^ fanh f^^"''"'^'"'^') cosh 

2LC{x',y>;a) ^C{x' ,y';aY/^ \ \ 2Ai / V 2Ai 

^ TT TT TT TT 

^3 = T f dx [ dy [ dx' [ dy' f25^^^(l - cosx)(l - cosx') 



4 C(27r)^ 



-TT — TT — TT — TT 



+52^^(1 - cos(x - y)){l - cos{x' - y')) S{x,y)H ( x,y,x',y'; 
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In G^{x, y), S^{x, y) and F there are additional corrections that come from inserting -ip into Gn{x, y) and S{x, y and 
Fq] then expanding these terms out to leading order in i]:' . For the purposes of illustration we will show 

G^{x,y) = kBTGnix,y) + {keTf [Af + [(2 - 2cos(a;) + (2 - 2 cos(y)]V'(ai sin^^o - 8a2 sin 2V'o)/(fcsr)] G^{x,yf. 

(E26) 



Appendix F. Derivation of self consistent approximation for periodic boundary conditions. 

Now, for the moment, let us concentrate on periodic boundary conditions and set r(x, y) = 0. We may represent 
each term in the perturbation theory diagrammatically. The Feynman rules; the correspondence between a graph and 
its algebraic expression are now given. On expanding out we need three types of diagrams (or vertices) to 

represent each 0"-term in this expansion, for n > 0. The Feynman rules are similar to those given in Appendix A of 
po|. but with a few modifications [s^. They are as follows. Each graph will contain Np vertices all of which will be 
connected to each other by lines. We assign a label i = 1 . . . Ny to each vertex. For each vertex i, representing (^", 
we must write down the following for a type 1 or 2 vertex 

(ai(-l)"/2-icos(V) + a2(-4)"/2cos(2^)) 5^ki^ki^+ki^_,+...ki/n\ when n is even 
(^ai(-l)("+i)/2-isin(^) +a2^^2 sin(2V.)j <5_,.,,.^+fe,_^+...fe. /n! when n is odd 



(Fl) 



and either one of the "form" factors 



n 

no 



l-e-^^™ , (F2) 



m—l 



depending on whether the vertex is type 1 or 2 , respectively. For a type 3 vertex, n is allowed only to be even, and 
we write down 

(ai(-l)"/2-i + 52(-4)"/2) <5-..,fc.„+^„_^+...4/n!, (F3) 

multiplied by the form factor 



no 



l_e-^(^™-!'™)). (F4) 



Each full Feynman graph will also consist of Ne external lines (connected to only one vertex) each associated with 
a wave vector qj = {xj /rQ,yj/ro,2Trnj/L) {j = 1, . . . ,Ne)- For each of these external lines we write down Gni(xi,yi) 
and set fc^ = qj for one of the fc^ — {x^^jr^, yl^/ r^, 27^711^ / L) in the vertex to which the line is connected. There will 
also be TV/ internal lines, each associated with a wave vector pk{k = 1, . . . , Nj), where each end is connected to two 
vertices, i and i' . For each of these internal lines we write down G{pk) and set fc^ = = pj for one of the k^^^ in 
each of the two vertices. Then all the wave vectors for the internal lines are summed over. Last of all, there is also 
a symmetry factor that multiplies this, which accounts for how many ways a term (graph) in the expansion may be 
generated. If we restrict the wave vectors to the form qj ~ {xj/r^, yj/ro, 0) we get back the results for the rigid body 
case ||2a| . 

To obtain the Hartree approximation we first consider the same set of graphs as for Appendix A of l20j as contri- 
butions to the full correlation function. The sum of these graphs we denote by Gin{x, y). These form a series which 
we may easily sum 

ksTGi^jx^y) _f27:ny f f L_ Ai ^ ^ 



(J \L J ^ V ^'^ vXi"'"' A~j j ((l-cos2;) + (l-cosy)) 

+ (a2.i (l-cos(x-2/)), (F5) 



26 



where 



.^i' ' \J ^ 2 (^aiCos(V') 



C 



„„„ I Aixl(£/(2Ai),a) 



4a2 cos(2V')exp (^-^^ 



(L/(2Ai),a) 
A„ 



C 



A9 1 I — (X — 



X2(-L/(2Ai),a) 



(F6) 



In the same procedure as discussed in Appendix A of ^(J|, we may replace Gn in each loop (of our diagrammatic 
series) with Gin- Now. Gi„ will be replaced on the l.h.s. of (F5) with a new correlation function 6*2,1. On the r.h.s of 
(F5) we replace Ai,i with Ai^2 and A2,i with A2,2- Expressions for Ai^2 and A2,2 are similar to (F6), but with A replaced 
by Ai_i and a replaced by ai = (Ai_i/A2,i)^. We then keep iterating this process until we have = Ai^oo = Ai_oo-ij 

^2 = X2.00 = A2,oo-i and uh = (Af^/A|^)^ so obtain Eq. (4.3) of the text. Through same reasoning as was discussed 
in Appendix. B of |2Qj, we find that ip = ipn where 



cos IpH 



ai 

4a2 V 



/3AiXi (i/(2Ai),, 



4A„ 



(F7) 



To calculate the free energy we consider the same set of graphs as those considered for free energy in [20j .To get 
the free energy in the Hartree approximation we then renormalize the sum of these graphs, by taking care in replacing 
Ai with \i , A2 with A|^, and a with an] and so arrive at Eq. (4.2) of the text. 



Appendix G. Correction to self consistent approximation for periodic boundary conditions from freely 

fluctuating ends. 

To calculate the corrections arising from free fluctuating ends we must utilize the results of both the previous 
two appendices. We perform a renormalization where we replace Ai, A2, a and V'o with Af^, Af'^, an and -0// in 
the results of Appendix E. The renormalization process relies on counter-terms which take care of terms already 
included in A{^, A^, an and 'i/'i?, to prevent over-counting. In the interests of brevity we will not discuss this process 
of renormalization, instead we refer those not acquainted with such processes to a standard text in field theoretical 
methods. 

The first quantity to consider is "0' which on renormalization becomes 

,/ ^ (ai sin-0g ~ 8a2sin20H) / /g\ 

(aicos?/'// - 4a2Cos2VjH) V ^7' ^ ' 



where 



Y" = / dx' I dy'il - cos(y')) 

f 4(27r)2ApL ' ' ^ ^ " 



L 



C(x',v'-aH) 2Ai(7(a;',y';aff)i/2 
^sinh cosh ^^(--^^^W)V-\ \ 



(G2) 



It is important to notice that there is no Fp^, this is because such a term is already accounted for in ipH- On 
renormalization, a counter term removes this term. 

Let us next consider G^(x, y), to obtain the correction from free boundary conditions we perform a renormalization 
of (E8) where we obtain 

G^^«(x, y) = keTG^ (x, y) + {kBT^G^ {x, yf 
(Af -f [(2 - 2cos(a;) + (2 - 2 cos(2/)]^/'^(ai sin^/>jj - 802 sin 27/-// )/(A;bT)) , (G3) 
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where 



Af = [(2 - 2 cos(a;)) + (2 - 2 cos(y))] A)'^ + (2 - 2 cos(a; - y))Aj^" , (G4) 



L/2 L/2 . 



2i C{x',v';aH) 



L/2 -L/2 

1^ 



4(7(:c',y';aH)i/2 II 2Af i I 2Af 



2i 



L/2 -L/2 



4C{x',y';aHy/^ \ \ 2Af \ 2Af 



and 



-2 (j = 1^— j +(Af) ((l-cosar) + (l-cosy)) 

+ (Af)"'(l-cos(x-y)). (G5) 

There is no A^, this is because such a term is already accounted for in G^{x,y), so is removed on renormahzation. 
We may go further and make the approximation that 

G^'«(x, y) ~ fcsTGf (x, y) + {ksTfC'^ {x, yf (Af + [(2 - 2 cos(a;) + (2 - 2 cos(y)]V'^ 
(ai sin V'ff - 8a2 sin 2V'ff)/(fcBT)) + {kBTfG^ (x, y)^ (Af + [(2 - 2 cos(x) + (2 - 2 cos(y))] 

V'ff(ai sinV-ff - 8a2sin2'!/'ff)/(fcBr)) + . . . 
~ kBTG^{x,y) + [kBTf (Af + [(2 - 2cos(x) + (2 - 2cos(y)]VH(ai sin^H - 802 sin2V'jf)/(A;i3T)) 

G^{x,y)Gl'^{x,y). (G6) 

Prom (G6) we may write 

kBTG^'^{x,y)-^ /27rnV ,,fx-2,,, , 

^ "^^ ^ = (^— j + (Af ) ((1 - cos x) + {l- cos y)) 

+ {XF)-\l-coB{x-y)), (G7) 

where 

(Af)"' = (Af)"'-2C-i (A;BrA)'^ + V4(aisinV'ff-8a2sin2V'jf)) 

(Af)"' = (Af)-' - 2C-'kBTA)'". (G8) 

Now, we renormahze S^{x,y) 

S'.''{..y) = S'-i^.y) + S-{.,yf (Sf - (°. - Sym2fe)feAf - + (I - cos,)] 



28 



-TT — TT 



' j dy' (^g[^\{2 - 2cos(x)) + (2 - 2cos(y)))(l - cos(a;')) 

— TT 

+5^^^(2-2cos(a;-y))(l-cos(a;'-y'))) (^a;, y, x', y'; aj/^ . 



There is no because such a term is already accounted for in S^{x,y). 

We may also renormalize the free energy, the free energy then becomes Ft = Fh + Fp + AF where 

F. = Jl^/h.ldyJ^{ ^[^^y^^-A ^'^ eoth { ^^-^y^-f'^^ 
2(277)2 y J ' I Af ^ C{x,y;l) ^ \ 

AF = -L{ilj'Hfmi - ksTG^, (GIO) 

and 



TT TT TT 



= \^^ j '^'^ j j '^^'Z '^y' [^9i\^-cosx){l-cosx') 

— TT — TT — TT — TT 

+52''^(l-cos(a;-y))(l-cos(a;'- J/')) S{x,y)H (^x,y,x' ,y'; -^,aH^ ■ (Gil) 

We may also calculate correlation functions within this approximation scheme. Here, we may use the results of the 
previous section, but we now make the following replacements: Ai is replaced by Af , a is replaced by ap = {^1/^2)^ 
and S{x,y) is replaced by S^'^{x,y). 

Appendix H. High Temperature Expansion. 

We start the analysis of this appendix by dividing (3.1) into three terms 

E[cp]=E^^, + E!^S]+ET[<f>], (HI) 



where 



Ei, = 3ao{Ri), (H2) 

L/2 

ELd<t>] = I rf^E E [(-ir«m(i?i) (cos (m (</.,, (^) - <^i+u-iW)) 

-L/2 ™=1 

+ cos(m {(pjiiz) - 0j_ii(z))) + cos(m (0j;(z) - (p.ji-i{z)))] , 

L/2 

'd(j)ji{z) 

-L/2 



ET[<t>]= j dzJ2^ 

T /O .1,1 



dz 



First, let us look at the high temperature expansion of the free energy. We start by expanding out the partition 
function in the following way 

- - (-0) n/-.<=. (- 4 ■ -) - (-1*) ■ <-) 
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Prom this expansion the lowest order terms in Z that depend on the coefficients a„ are 



ym m 

1 ~ 2{kBT)^ 



exp - 



-^int 



L/2 L/2 



) n / ^^"i(') I '^^ii I ^4 EE J J 

i,l f ,1' .7'",'" r /o r /o 



-L/2 -L/2 



>s (rn. - rj,-n'i^) + | (7,'/' - 7i'-U') + 4^' - " 

cos (m (^fj,i,{z) - 4>],i,_^{z) + J {'yyv - Ij'V-i) + - j'-i)) " 
cos (to (4>^j,i,iz) - 4>f+ii>-iiz) + Y hj'i' - 7j'+K'-i) + 4>fi' - (pf+w-i)) 



cos ( m ( ^„i„{z') - ^„_^i„{z') + Y ilyi" - 7j"-U") + 0°"; 



+ COS ( TO ( ^„i„{z') - ^„i„_-i_{z') + — - + - 



cos ( TO ( 4fy,i„{z') - 4fj„^-yi„_-^{z') + y (7j"i" - 7j/, + i;„_i) + </)°„;„ - (j)P„ _^_^i„ 



exp 



(H4) 



where we have employed (3.7) of the text but, now separating out the rigid body mode so that (j)^^ (z) = (j)^^ {z) + <fPy , 

and the spatial average of 4fji{z) along the length of the molecule is zero. We find that the only a few of these terms 
survive integration over the rigid body modes. 



m m 



L/2 L/2 



n/p^,^,(.)/d7,.x: J dz J d.'exp(-^) 



-L/2 -L/2 



cos 



(m (4,(z) - - + rr_,,{z') + {z- z') {l,v - li'-w)/L)) + 

cos (to {f.,,{z) - ^,,{z') - ^,,_,{z) + ^,,_,{z') +{z- z') i^jn' - Jjn'-i) /L)) 
+ COS (to (^,,{z) - ^,,{z') - ^,+i,_i(z) + ^,+n'-ii^') Hz - A iJj'V - ir+ii'-i) im . 
Now we may rewrite Et [0] as Et = E^ + E^ where 



3,1 3,1 \ / 



We find that we may then write the following expansion of In Z 

lnZ~lnZo + ^Zr/Zo 



(H5) 



(H6) 



(H7) 



L/2 L/2 



Zo 4(fcs 



3 _i/2 _L/2 



X 

+ [cos 
+ (cos 



s (to U^.,,{z) - ^,,{z') - 4_,,(z) + (cos(to((z - z') ijjn' - l3'-w)IL))) 

(to (4,(z) - - ^'v-x{z) (cos (TO {{z - z') {^rv - ^3'i'-i) / L))) e', 

(to (^^,,(z) - ~4>^y,{z') - 4+i,_i(z) + 4+i,_i(^'))))^, (cos (TO {{z - z') {^yi, - ^y+n'-l) /L)))^^ 



30 



where 



(m (4,(z) - - + 

n cos (m - ^,,{z') - ^,_,,{z) + exp (-^) 

J,' 

... n/rf7..cos(.n((f^(,,.-,,,_,0))exp(-S^) 
cos m (ciyv - 7,v_ip) \\\ = 



(H8) 



The term (cos (n {{z — z') (7^'/' — 7j'_i;') / L))) is easy to evaluate and we find 



cos ( m ( If_^ (^^,,, _ 7j/^ir)) n = (cos ( m ( ^^-^ (7j'/' - 7j'/'-i) 



cos m 



The term {^cos — 0^,;,(z') — 0^, + 0^, ) p slightly more involved, where to evaluate this 

term we express (z) in terms of its Fourier series (j)^i{z) ~ J2 ^ji '^'^P ( ^^'^"^ ) we then find that 

,(z) + iP, .„(z')]]) 

E\ 

The sum may be evaluated leaving us with 

(cos (m (f^,,{z) - ^„{z') - ^,^,,^,{z) + 0^,+,,„,(z))))^^ = exp (-^^^^^^) exp • (Hll) 

Putting this all together we find 

L/2 L/2 

3a^iV f f ( \z^ z'W _ial^L'^N fm^L 



(m (0^,,(z) - - + 4_i,(z') 

(m (0^,,(z) - - 4+i,_i(z) + 



-L/2 -L/2 

where 

[8x + 8exp(-x/2)~8] 
/2(.a;j = 2 • (J^-Ldj 

If we truncate the sum in (HI) over the am, according to for a DNA molecule, at m = 2, then we may also 
calculate the next to leading order correction to the free energy per molecule, 

(TJX A 

(H14) 



3a2alL^h,2iL/Xp) 



8(fcsT)2 
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where Fq is a term independent of ai and 02, and 

24[exp(-x/2)(4 + x) + (-4 + x)] 
/3,iW = ^3 , 

^ , , [exp(-a;/2)(128 + 24a;) + 36a;-126-2exp(-2a;)] 

/3,2(a:j = ^-3 . (H15j 

We may look at limiting cases. When the molecules are either very short or rigid, so that L <C Ap, then we find 
that we may write 

^ = ^°-^(fc^T) 2(feT)^ +^77r^ + ---- (H16) 

This yields the following specific heat 

Here we have an additional term when compared with [20], this term was overlooked previously. We may also look 
at the limit where molecules are very long or soft. Here we find 

When comparing the two limits we see that as we move from short rigid molecules to long flexible molecules our 
high temperature expansion changes from a power series expansions in La\/{kBT) and La2/{kBT) to power series 
expansions in ApOi/ (ksT) and Xpa2/{kBT). So for L 3> Ap, the high temperature expansion being valid when we have 
Xpai/{kBT) <C 1 and Xpa2/{kBT) <C 1, not necessarily Lai/{kBT) <C 1 and La2/{kBT) <C 1. It is also interesting to 
look at the correlation function 

(exp (in iM^) - 'l^i'i'iz')))) = I n/ ^'^(^) «'^P ('^i'(^) - '/'/H^'))) exp (-^) (H19) 

evaluated for nearest and next nearest neighbors. To generate the high temperature expansion we make the following 
expansion 

f\ -1 




In general the leading order term in the correlation function- that describes how azimuthal correlations are lost- is of 
the form 



{e^^{in{M^)-4>Mm)^n, (^^) /, (z'i'^) ^^^^^ 

Here, b is the smallest number of links or bonds between sites and {i' A link being defined as a translation 
of ±i?iu, ±i?iv or ±i?i (v — u) between two sites on the 2-D lattice. The factor nt is the number of paths using the 
smallest number of links that go between the two sites and This is illustrated for the examples that we 

give explicit results in Fig. 8. Also the scaling function should always have the property that for torsionally rigid 
molecules /b = 1. 

First let us calculate the correlation function where 6=1. This is the correlation function between nearest neighbors. 
We find to leading order in the high temperature expansion. 

{exjp{in{(l)ji{z) - (j)j-u{z')))) = {exp{in{(l>ji{z) - 0^7-1(2')))) 
= (exp(m(<^,,(^) - cj>j+u-i{z'm = ^i^P/i (l' X= ^) . (H22) 
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/ % — m^: % . m 

si w*^ k 

FIG. 8: The paths with the shortest number of hnks for the correlation functions a.) (cxp (ire (<^j;(z) — b.) 
(exp(m(0j;(2:) - (?ij+i;(2;')))) c.) (cxp(m((/)jj(«) - (?!)j+2;(2')))) and d.) (cxp (in - 0j+ii+i(2:')))). Now, a.) and b.) 

contain only one link, while c.) and d.) two. d.) is the only one, shown here, to have two different paths, path I and II. 



where 



fi{x,x';y) = i (^4exp(^-^^^^^ + y\x - x'\exp 
-2exp (^--J I exp 



y\x — X 



4 

y{x + x') 



(H23) 



When L we do, indeed, find that fi{x,x';y) ~ 1. 



We may also compute correlations where 6 = 2. This is between next to nearest neighbors and next to next to 
nearest order correlations. We find to leading order 

2 {exp{in{(j)ji{z) - (j)j_2i{z')))) = 2 (exp {in{(l)ji{z) - (j)ji_2{z')))) 
= 2 (exp {in{(j>ji{z) - (j)j+2i-2{z')))) = {exp {in {(pji{z) ~ (pj_ii_i{z')))} 
= {exp {in {4>ji{z) - (t>j_2i+i{z')))) = {exp{in{<j)ji{z) - (j)j_ii+2{z')))) 



an{-iy 
2kBT 



where 



h {x,x';y) = 



4f{x,x';y) , 8 



1 1 
-+4- 



exp 



2(x + x') 
H — ^ exp 



(-!) 



exp 



exp 



y 

y{x + x'y 
4 

y{x + x') 



+ ^ exp - 



. z z' n L 



y\x- x'\ 



(H24) 



exp 
exp 



+ 



2\x- x'\ {x - x'f 



+ 



y 2 

y{x - x') 



exp 



4 

y{x + x') 

4 

y{x + x') 

4 

y\x — X' 

4 

y{x - x') 



Again, when \ ^ L, we find that f2{x,x';y) ~ 1. 

We may also look at self correlations (6 = 0), due to torsional fluctuations, for a single DNA molecule 

{exp {in {(p J I (z) - (j)ji{z')))) = exp 



(H25) 



-n^ \z — z'\ 



4A„ 



(H26) 
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